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We propose the Bose-Einstein condensation (BEC) and superfluidity of quasi-two-dimensional 
(2D) spatially indirect magnetobiexcitons in a slab of superlattice with alternating electron and 
hole layers consisting from the semiconducting quantum wells (QWs) and graphene superlattice 
in high magnetic field. The two different Hamiltonians of a dilute gas of magnetoexcitons with a 
dipole-dipole repulsion in superlattices consisting of both QWs and graphene layers in the limit of 
Ch ' high magnetic field have been reduced to one effective Hamiltonian a dilute gas of two-dimensional 

excitons without magnetic field. Moreover, for A'^ excitons we have reduced the problem of 2 A*' x 2 
dimensional space onto the problem of A'^ x 2 dimensional space by integrating over the coordinates 
' of the relative motion of an electron (e) and a hole (h). The instability of the ground state of 

the system of interacting two-dimensional indirect magnetoexcitons in a slab of superlattice with 
alternating electron and hole layers in high magnetic field is established. The stable system of indi- 
' ■ rect quasi-two-dimensional magnetobiexcitons, consisting of pair of indirect excitons with opposite 

(— I ' dipole moments is considered. The density of superfiuid component ns{T) and the temperature of 

I , the Kosterlitz-Thouless phase transition to the superfiuid state in the system of two-dimensional 

^ ■ indirect magnetobiexcitons, interacting as electrical quadrupoles, are obtained for both QW and 

■ ' graphene realizations. 

PACS numbers: 71.35.Ji, 71.35.Lk, 71.35.-y, 68.65.Cd, 73.21.Cd 
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I. INTRODUCTION 



B 

Q ■ The many-particle systems of the spatially-indirect excitons in coupled quantum wells (CQWs) in the presence or 
Q ■ absence of a magnetic field B have been the subject of recent experimental investigationsii^i^i^. These systems are 
of interest, in particular, in connection with the possibility of Bose-Einstein condensation (BEC) and superfluidity of 
indirect excitons or electron-hole pairs, which would manifest itself in the CQW as persistent electrical currents in 
each well and also through coherent optical properties and Josephson phenomenai^i^i^iiSii^ In high magnetic fields, two- 
I dimensional (2D) excitons survive in a substantially wider temperature range, as the exciton binding energies increase 
^\ , with magnetic field i-'^^'^^i^^'^'^i-'^'^'^^i^^ The problem of essential interest is also collective properties of magnetoexcitons 
' in high magnetic fields in superlattices and layered systemiii 

' In this Paper we propose a new physical realization of magnetoexcitonic BEC and superfluidity in in superlattices 

■ with alternating electronic and hole layers, that is in a sense representing array of CQWs or graphene layers with 
spatially separated electrons (e) and holes h with spatially separated electrons and holes in high magnetic field. 

00 : Recent technological advances have allowed the production of graphene, which is a 2D honeycomb lattice of carbon 

■ atoms that form the basic planar structure in graphit o^^'^^ Graphene has been attracting a great deal of experimental 
^ ' and theoretical attention because of unusual properties in its band structure22i2ii2^i2^. It is a gapless semiconductor 

with massless electrons and holes which have been described as Dirac-fermiona^l. Since there is no gap between the 
. ^ conduction and valence bands in graphene without magnetic field, the screening effects result in the absence of excitons 
^ . in graphene in the absence of magnetic field. A strong magnetic field produces a gap since the energy spectrum becomes 
■ ■ ' discrete formed by Landau levels. The gap reduces screening and leads to the formation of magnetoexcitons. We also 
consider magnetoexcitons in the superlattices with alternating electronic and hole graphene layers. We suppose that 
recombination times can be much greater than relaxation times due to small overlapping of spatially separation of 
e- and h- wave functions in CQW or graphene layers. In this case electrons and holes are characterized by different 
quasi-equilibrium chemical potentials. Then in the system of indirect excitons in superlattices, as in CQW^'-, the 
quasiequilibrium phases appear—. While coupled-well structures with spatially separated electrons and holes are 
typically considered to be under applied electrical field, which separates electrons and holes in different quantum 
wells^ii^, we assume there are no external fields applied to a slab of superlattice. If "electron" and "hole" quantum 
wells alternate, there are excitons with parallel dipole moments in one pair of wells, but dipole moments of excitons 
in another neighboring pairs of neighboring wells have opposite direction. This fact leads to essential distinction of 
properties of e — /i system in superlattices from ones for coupled quantum wells with spatially separated electrons and 
holes (where indirect exciton system is stable due to dipole-dipole repelling of all excitons). This difference manifests 
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itself already beginning from three-layer e — h — eovh — e — h system. We assume that alternating e — h — e layers 
can be formed by independent gating with the corresponding potentials which shift chemical potentials in neighboring 
layers up and down or by alternating doping (by donors and acceptors, respectively). 

In this Paper we reduce the problem of magnetoexcitons to the problem of excitons at i? = 0. The unstability 
of the ground state of the system of interacting indirect excitons in slab of superlattice with alternating e— and h— 
layers is established in strong magnetic field. Two-dimensional indirect magnetobiexcitons, consisting of the indirect 
magnetoexcitons with opposite dipole moments, are considered in high magnetic field. The radius and the binding 
energy of indirect magnetobiexciton are calculated. These magnetobiexcitons repel each other as electrical quadrupoles 
at long distances. In result, the system of indirect magnetobiexcitons is stable. In the ladder approximation collective 
spectrum of the weakly interacting by the quadrupole law two-dimensional indirect magnetobiexcitons is considered. 
The superfluid density ng{T) of interacting two-dimensional indirect magnetobiexcitons in superlattices is calculated 
at low temperatures T. We analyze the dependence of Kosterlitz-Thouless transition^^ temperature and superfluid 
density on magnetic field. 

The rest of the Paper is organized as follows. In Sec. HIl we derive an effective Hamiltonian of magnetoexcitons 
in both CQWs and two graphene layers in strong magnetic field. In Sec. IIIH we prove the instability of dipole 
magnetoexcitons in QW and graphene superlattices due to the attraction of oppositely directed dipoles. BEC and 
superfluidity of quadrupole magnetobiexcitons in QW and graphene superlattices has been established and analyzed 
in Sec. IIVI We present and discuss the numerical results in Sec. |Vl 



II. EFFECTIVE HAMILTONIAN OF MAGNETOEXCITONS IN STRONG MAGNETIC FIELD 



We start with two interacting excitons in CQW in the presence of the external magnetic field B. Without loss of 
generality we can take the magnetic field to be in the z direction. We find it convenient to work with the symmetric 
gauge for a vector potential Ae(/j) = 1/2 [B x re(/i)]. We employ Jacobi coordinates for the exciton center of mass 
R = [rUpVe + mhXh) / i'm-e + 'm-h) and relative motion of electrons (e) and holes (h) r = re — and in these coordinates 
the Hamiltonian H for 2D spatially separated excitons is 



iJ = dR. dr 




^JdRiJ dn J dK2 J dr2V'^(Ri,ri)?/;t(R2,r2) 

i^j—e^h 



(1) 



where {R, r) and ip(R, r) are the creation and annihilation operators for magnetoexcitons; and are electron and 
hole locations along quantum wells, respectively; D is the distance between electron and hole quantum wells, e is the 
charge of an electron; c is the speed of light and e is a dielectric constant (for graphene layers e is the dielectric constant 
of the matrix in which graphene layers are embedded). In Eq. ^ we use the two-particle potentials for the electron- 
electron, hole-hole, electron-hole and hole-electron interaction, for electrons and holes from two different excitons: 
U'^rei - r,2) = eV(e|re 
f/''"(rM - v,2) = -eV(ev4^ 



re2|), [/''''(rw - Vh2) = eV(e|rM - Vh2\), U'H^ei - rh2) = -e^ I - Vh2? + D"^ 



■V,2?+D^). 

The four-component Hamiltonian of an isolated electron-hole pair in bilayer graphene with spatially separated 
electrons (e) and holes [h) in one valley in magnetic field B neglecting the Coulomb interaction is given by22 



H = vp 



I Px^ +ipy^ 

pt^-ip-y^ 

pl'^' - ^p'^' 
(h) . ■ (h) 

Px +iPy 







(2) 



where 



p(-) = -ihVe + -Ae; 

c 



(3) 
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where e is an electron charge; c is the speed of hght; and are two-dimensional (2D) vectors of coordinates of 
an electron and hole, correspondingly; Ae and A/j are the vector potential of an electron and hole, correspondingly; 
vp = ^/2>at/{2h) is the Fermi velocity of electrons in graphene (a = 2.566A is the lattice constant, t fa 2.71eV is the 
overlap integral between the nearest carbon atoms }^ . 

Note that the ratio of the contribution to the energy from the Zeeman term AEz{B) to the characteristic separation 
between the nearest Landau levels AEl{B) is negligible (for quantum wells: this ratio is AEz{B)/ AEl{B) w 



(0) 



0.1 ; for graphene layers: at B ^ lOT this ratio is AEz{B) / AEl{B) w ^bB / {V^hvprB^) w 5 x 10" 



where /is = he/{2m^e^c) is the Bohr magneton; 

,(0) 



0.1mi°^ is the effective electron mass in GaAs/AlGaAs 



quantum wells; rUe is the mass of a free electron). Therefore, the contributions to the single-electron Hamiltonian 
from the Zeeman splitting set identically to zero analogously to Refs. j23ll27[ . We take into account the energy 
degeneracy corresponding to two possible spin projections in quantum wells and graphene and two graphene valleys 
(two pseudospins). Since electrons on a graphene lattice can be in two valleys, there are four types of excitons in 
bilayer graphene. Due to the fact that all these types of excitons have identical envelope wave functions and energies^!, 
we consider below only excitons in one valley. Also, we use no = n/(4s) as the density of excitons in graphene 
superlattice, with n denoting the total density of excitons, s is the spin degeneracy (equal to 4 for magnetoexcitons 
in bilayer graphene). Besides, we use tiq — n/ s as the density of excitons in quantum wells, with n denoting the total 
density of excitons, s is the spin degeneracy (equal to 4 for magnetoexcitons in coupled quantum wells). 

A conserved quantity for an isolated electron-hole pair in magnetic field B is the exciton generalized momentum P 
defined as 

P = -ifNe - iK^h + - (Ae - Ah) - -[B X (re - r^)] 

c c 

for the Dirac equation in graphene layers^! as well as for the Schrodinger equation in CQW a"'^°'^^i^^ . 

The Hamiltonian of a single isolated magnetoexciton without any random field (ye(re) — Vh{vh) = 0) is commutatcd 
with P, and hence they have the same the eigenfunctions, which have the following form (see Refs. [Tolf29j V 



*fcp(R,r) 



expjzR (p + |BxR) +z7^|$,(P,r), 



(4) 



where $fc(P,r) is a function of internal coordinates r and the eigenvalue P of the generalized momentum, and k 
represents the quantum numbers of exciton internal motion. The wave function of the relative coordinate for e and 
h spatially separated in different graphene layers <l>(r) can be expressed in terms of the two-dimensional harmonic 
oscillator eigenfunctions $„j_„2(i")- I^or an electron in Landau level and a hole in level n_, the four-component 
wave functions for the relative coordinate are^ 



.(r)-(V2)' 



V 



-*|«+|-i>-|(i-) 
-*|n+|,|«_|-i(r) 
*|™+|,|n_|(r) 



L(r)\ 



(5) 



where s± = sgn(n±). The corresponding energy of the electron-hole pair Enl,n- (which is the eigenvalue of the 
Hamiltonian (Eq. ([2]))) is given by^^ 



<\n_ = — V2 sgn(n+)yM -sgn(n-)v^ 



(6) 



where rs = ^/ ch/{eB) is a magnetic length; vp — y/3at/{2fi) is the Fermi velocity of electrons in graphene (a = 2.566A 
is a lattice constant, t w 2.71eV is the overlap integral between the nearest carbon atoms).— The two-dimensional 
harmonic oscillator wave functions eigenfunctions corresponding to e and h in spatially separated CQWs 4>„j_„2(r) 
are given by^l 



,(r) = (27r)-i/22-l™l/- 



. hi 1 . .^rl"! 

Jd 



exp 



-imd) s- 



L, 



\m\ 



2r2 



(7) 
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where denotes Laguerre polynomials; m = n-i — 77.2; n = minfrti , na), and sgn(TO)'" — s- 1 for m = 0. The 

expression for <I>fc(P, r) for CQW is given in Ref. [Tolf3 |. and in Ref. [13]) for graphene layers. In high magnetic fields 
the magnetoexcitonic quantum numbers k = {n_|_,n_} for an electron in Landau level and a hole in level n_; 
7 = (m/i - me)/(TO?i + me). 

In a strong magnetic field at low densities { n <^ r^^ ) (rs = ^/hc/ ieB) is the magnetic length) indirect magne- 
toexcitons repel as parallel dipoles, and we have for the pair interaction potential: 

t/(|Ri - R2I) = U'' + ii'"'' + + U''' ~ n - 

e|Ki — K,2| 

If we expand the magnetoexciton field operators in a single magnetoexciton basis set \I'fcp(R, r): ■0^(R, r) — 

J2kP ^fcp(-^' '')^Ip' ^(R-jI") — Sfcp ^fcpCR-i r)afcp, where a[.p and a^p are the corresponding creation and annihila- 
tion operators of a magnetoexciton in (fc, P) space and substitute the expansions for the field creation and annihilation 
operators into the Eq. ([T]) and obtain the effective Hamiltonian in terms of creation and annihilation operators in P 
space. In high magnetic field, one can ignore transitions between Landau levels and consider only the lowest Landau 
level states (for CQWs n+ = n_ = 0; for graphene layers n+ = n_ = 1). 



Using the orthonormality of the functions $fe(0,r) we obtain the effective Hamiltonian HcS in strong magnetic 
fields. 

Due to the orthonormality of the wave functions <I>„_|_^„_ (0, r) the projection of the Hamiltonian (Eq. ([1])) onto 
the lowest Landau level results in the eff'ective Hamiltonian, which does not refiect the spinor nature of the four- 
component magnetoexcitonic wave functions in graphene. Since typically, the value of r is rs, and P Tt-I^b in this 
approximation, the effective Hamiltonian i/eft in the magnetic momentum representation P in the subspace the lowest 
Landau level (for QWs = n_ = 0; for graphene layers n+ = n_ = 1) has the same form (compare with Ref.jll) as 
for two-dimensional boson system without a magnetic field, but with the magnetoexciton magnetic mass jub (which 
depends on B and Z); see below) instead of the exciton mass {M = me -t-m^), magnetic momenta instead of ordinary 
momenta and renormalized random field (for the lowest Landau level we denote the spectrum of the single exciton 
eo(P) = eoo(P)): 



ffcff = ^ £o(^')apap + ^ 



p 

X 



^ (Pi,P2 P3,P4)4iakaP3«P4, (8) 

Pl,P2,P3,P4 



where the matrix element ^Pi,P2 C/|P3,P4y is the Fourier transform of the pair interaction potential U{R) — 

e^D^ /tB?. For an isolated magnetoexciton on the lowest Landau level at the small magnetic momenta under con- 
sideration, eo(P) ~ -^^/(Swb), where is the eff'ective magnetic mass of a magnetoexciton in the lowest Landau 
level and is a function of the distance D between e - and h - layers and magnetic field B (see Ref.fl^). In strong 
magnetic fields at D ^ the exciton magnetic mass is mB{D) = ^D'^ / {^'^f%) for QWs^ and mB{D) = eD^ / (ie^rg) 
for graphene layers^. 

Note that the ratio of the contribution to the energy from the Zeeman term AEz{B) to the characteristic separation 
between the nearest Landau levels AEl{B) is negligible (for quantum wells: this ratio is AEz{B) / AEl{B) « 

melra-e^ « 0.1 ; for graphene layers: at B = lOT this ratio is AEz{B) / AEl{B) « ^iBB/{V2hvFrB^) ~ 5x IQ-^, 

where /is = he/ (2m^^c) is the Bohr magneton; irie ~ O.lmi^'' is the effective electron mass in GaAs/AlGaAs quantum 

wells; mi°'' is the mass of a free electron). Therefore, the contributions to the single-electron Hamiltonian from the 
Zeeman splitting and very small pseudospin splitting in graphene (caused by two valleys in graphene) set identically 
to zero analogously to Refs. p3ll27j . We assume the energy degeneracy respect to two possible spin projections in 
quantum wells and graphene and two graphene valleys (two pseudospins) . Since electrons on a graphene lattice can 
be in two valleys, there are four types of excitons in bilayer graphene. Due to the fact that all these types of excitons 
have identical envelope wave functions and energies^I, we consider below only excitons in one valley. Also, we use 
no = n/(4s) as the density of excitons in graphene superlattice, with n denoting the total density of excitons, s is 
the spin degeneracy (equal to 4 for magnetoexcitons in bilayer graphene). Besides, we use hq = n/s a.s the density 
of excitons in quantum wells, with n denoting the total density of exciton, s is the spin degeneracy (equal to 4 for 
magnetoexcitons in coupled quantum wells). 

For large electron- hole separation D ^ vb, transitions between Landau levels due to the Coulomb electron-hole 
attraction can be neglected, if the following condition is valid, (for QWs Ei, = e^ /{tbD) ^ hujc — heB{me + 
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irih) / {2meTnhc); for graphene layers Ef, = 4e^/(eD) <^ Uvf/tb', where Ef, and Wc are the magnetoexcitonic binding 
energy and the cyclotron frequency, correspondingly). This corresponds to high magnetic field B, large interlayer 
separation D and large dielectric constant of the insulator layer between the graphene layers. In this notation, 
vf — y/iat/[2K) is the Fermi velocity of electrons. Also, a — 2.566A is a lattice constant, t « 2.71eV is the overlap 
integral between nearest carbon atoms^^. 



III. INSTABILITY OF DIPOLE MAGNETOEXCITONS IN QW AND GRAPHENE SUPERLATTICES 

Let us show the low-density system of weakly interacting two-dimensional indirect magnetoexcitons in superlattices 
is instable, contrary to two-layer system in CQW. At small densities nrg <C 1 the system of indirect excitons at low 
temperatures is the two-dimensional weakly nonideal Bose gas with normal to wells dipole moments d in the ground 
state {d ~ eD, D is the interwell separation), increasing with the distance between wells D. In contrast to ordinary 
excitons, for low-density spatially indirect magnetoexciton system the main contribution to the energy is originated 
from dipole-dipole interactions C/_ and J7+ of magnetoexcitons with opposite and parallel dipoles, respectively. Two 
parallel (-1-) and opposite (— ) dipoles in low-density system interact as U+ = —U- = e^D^/eR^, where e is the 
dielectrical constant; R is the distance between dipoles along wells planes; we suppose that D/R <C 1 and L/R 1 
(L is the mean distance between dipoles normal to the wells). We consider the case, when the number of quantum 
wells k in superlattice is restricted k ^ 75^7^1 and this is valid for small k or for sufhciently low exciton density (n 
is exciton surface density). 

The distinction between magnetoexcitons and bosons manifests itself in exchange effects. ^•■^^ The exchange inter- 
action in spatially separated system is suppressed in contrast to e — ft. system in one well due to smallness of tunnel 
exponent T connected with the penetration through barrier of dipole-dipole interaction. Hence, at 3> rB exchange 
phenomena, connected with the distinction between excitons and bosons, can be neglected for the both QWs and 
graphene layers"^°. 

For the analysis of stability of the ground state of the weakly nonideal Bose gas of indirect excitons in superlattices 
we apply the Bogolubov approximation. The total Hamiltonian Htot of the low-density system of indirect excitons 
in superlattice is given hy.Htot — Hq + Hint- Here Hq is the effective Hamiltonian of the system of noninteracting 
magnetoexcitons: 

Ho = Yl ^o(p)(a+ap + 6+5p + atpa_p + b+pb_p), (9) 
p 

where £o(p) = P^/i'^mB) is the spectrum of isolated two-dimensional indirect magnetoexciton; p represents the 
excitonic magnetic momentum, a^, Up, b^, bp are creation and annihilation operators of magnetoexcitons with up 

and down dipoles; Hint is the effective Hamiltonian of the interaction between magnetoexcitons: 

Hint = ^ I]pi+p2=p3+P4("p4"p3'^P2«Pi + 
^lt4^P3^P2Vl - «p4«lt3V2V - ^P4^P3«P2«P1- (10) 

S is the surface of the system. Let us consider the temperature T = 0. Assuming the majority of particles are in 
the condensate {N — Nq/Nq <C 1, where N and Nq are the total number of particles and the number of particles in 
condensate), we account as in Bogolubov approximation only interaction between condensate particles and excited 
particles with condensate particles, neglecting by the interaction between noncondensate particles. Then the total 
Hamiltonian transforms to: 

Htot = ^Y.p^o[>^o{p){a+ap + b+bp + atpa_p+ 

b-pb-p) - Un{a+btp + apb_p + atpb+ + a_pbp (11) 
+4^P + a-pb-p + apb+ + a_pbtp)], 

In Eq. (|lip terms, arising from first and second terms of the Hamiltonian Eq. (|10p . which describe the repulsion of the 
indirect magnetoexcitons with parallel dipole moments, are compensating by other terms of the Hamiltonian Eg. ljlOp . 
describing the attraction of indirect magnetoexcitons with opposite dipoles. In the result only terms describing the 
attraction survive. Let us diagonalize Hamiltonian Htot by using of the unitary transformation of the Bogolubov 
type^ 



+ 



«p = -/T=^=^(^p + ^pO'-p + BpP^p, 
bp = ^i_il_gl (/3p + ApPtp + Bpatp) 



(12) 
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where the coefRcients Ap and Bp are found from the condition of vanishing of coefficients at nondiagonal terms in 
Hamiltonian. In result we obtain 

Htot = £(P)("P«P + f^pf^p) (13) 

p^O 

with the spectrum of quasiparticles e{p): 

e{p) = VieoipW - {nUf- (14) 

At small momenta p < y/2mBnU the spectrum of excitations becomes imaginary. Hence, the system of weakly 
interacting indirect magnetoexcitons in slab of superlattice is unstable. It can be seen that the condition of the 
instability of magnetoexcitons as stronger as magnetic field higher, because niB increases with the increase of magnetic 
field, and, therefore, the region of p resulting in the imaginary collective spectrum increases as B increases. 

This, on first view, strange result can be illustrated by the following example. There are equal number of dipoles 
oriented up and down. Let us consider four dipoles, two of them being oriented up and two — down. It is easy to 
count that number of repelling pairs is smaller than that of attracting ones. The prevailing of attraction leads to 
unstability. 



IV. BEC AND SUPERFLUIDITY OF QUADRUPOLE MAGNETOBIEXCITONS IN QW AND 

GRAPHENE SUPERLATTICES 



Let us consider as the ground state of the system the low-density weakly nonideal gas of two-dimensional indirect 
magnetobiexcitons, created by indirect magnetoexcitons with opposite dipoles in neighboring pairs of wells (Fig. [1]). 
The small parameter at the adiabatic approximation is the numerical small parameter which is equal to the ratio of 
magnetobiexciton and magnetoexciton energies or the ratio between radii of magnetoexciton and magnetobiexciton 
along quantum wells (or graphene layers) (see, e.g., Ref. [SS*]). These parameters are small, and they are even smaller 
than analogous parameters for atoms and molecules. The smallness of these parameters will be verified below by 
the results of the calculation of indirect magnetobiexciton. Here it was assumed, that the distance between wells (or 
graphene layers) D is greater than the radius of indirect magnetobiexciton ai,. The potential energy of interaction 
between indirect magnetoexcitons with opposite dipoles U{r) has the form shown on Fig. [2](r is the distance between 
indirect magnetoexcitons along quantum wells/graphene layers): 



C/ r = + , (15) 

er eVr^ + eVr^ + ^D^ 

At r > l.llD indirect magnetoexcitons attract, and at r < l.llD they repel. The minimum of potential energy U{r) 
locates at r = ro w 1.67D between indirect excitons. At large D one can expand the potential energy U{r) on the 
parameter r — ro/D <^ 1: 

2 2 

C/(r) = -0.04^ +0.444^(r-ro)2. (16) 

So at large D magnetobiexciton levels correspond to the two-dimensional harmonic oscillator with the frequency 
w = 0.88eV(mBeL'3): 

S„ = -0.04^ + 2x/2So f^j (" + 1), (17) 

where Eq = tubc'^ /{ti^e), r* — h'^e/{2mBe'^), e'^ = 0.88e^. In the ground state the characteristic spread of magneto- 
biexciton Ofc along quantum wells/graphene layers (near the mean radius of magnetobiexciton vq along wells/graphene 
layers) is: 



where a^^ = (Sr'^^y^^D^^^; r'^^ = h^e/{2mBe^) is the two-dimensional effective Bohr radius with the effective magnetic 
mass ruB- Hence, the ratio of the binding energies of magnetobiexciton and magnetoexciton is: Ebex/Egx = 0.04 <ti 1 
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FIG. 1: Two-dimensional indirect magnetobiexcitons consisting of indirect magnetoexcitons with opposite dipole moments, 
located in neighboring pairs of wells/graphene layers. 
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FIG. 2: The potential energy U{r) of the interaction of indirect magnetoexcitons with opposite dipoles, located in neighboring 
pairs of wells or graphene layers (in units of the binding energy of indirect magnetoexciton Eex ~ e^/eD), as a function of the 
distances r between magnetoexcitons along the wells/graphene layers (in units of D). 



at D ^ ttex (the ratio of radii of magnetoexciton and magnetobiexciton is rQ/a^x = 0.67(8r^^)^/'*Z?^^/'* ^ 1). So 
adiabatic condition is valid. 

The mean dipole moment of indirect magnetobiexciton is equal to zero. However, the quadrupole moment is nonzero 
and equal to Q = 3eD^ (the large axis of the quadrupole is normal to quantum wells/graphene layers). So indirect 
magnetobiexcitons interact at long distances R:^ D as parallel quadrupoles: U{R) — 9e^D^ / (eR^). 

Exchange effects, connected with the distinction between low-density indirect magnetobiexcitons and bosons, can 
be suppressed due to the negligible overlapping of wave functions of two magnetobiexcitons on account of the po- 
tential barrier, associated with the quadrupole repulsion of indirect magnetobiexcitons at long distances analogously 
to dipole magnetoexcitons. At large D the small tunnelling parameter connected with this barrier has the form 
exp[— 0.93 y^UJr^]. Hence, at _D 3> r*^^ exchange effects for indirect magnetobiexcitons can be neglected. 

We account the scattering of magnetobiexciton on magnetobiexciton by using of the results of the theory of two- 
dimensional Bose-gas;^ The chemical potential /i of two-dimensional biexcitons, repulsed by the quadrupole law, in 
the ladder approximation, has the form (compare to Refs. iSflSS]): 



m% log ?i'*/\2/3/ (8TT{9m%e^D'^)^/^nbe 



(19) 



where Ubex = n/{2s) is the density of magnetobiexcitons in QWs and ribex ~ n/{8s) in graphene layers; = 2mB 
is the mass of a magnetobiexciton. 

At small momenta the collective spectrum of magnetobiexciton system is the sound-like e{p) — CsP {cg — \J~^iJm}^ 

is the sound velocity) and satisfied to Landau criterion for superfluidity. The density of the superfluid component 
ns{T) for two-dimensional system with the sound spectrum can be estimated as)^ 



(20) 
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FIG. 3: Dependence of temperature of Kosterlitz-Thouless transition Tc = Tc(B) (in units K) for tlie superlattice consisting 
of quantum wells (for GaAs/AlGaAs: t = 13); and (for graphene layers separated by the layer of Si02'- e = 4.5) on the 
magnetoexciton density n (in units cm~^) at the interlayer separation D = Wnm at different magnetic fields B: B = 20T - 
solid curve for quantum wells and dotted curve for graphene layers; B = 15T - dashed curve for quantum wells and dashed- 
dotted curve for graphene layers; B = IQT - thin solid curve for quantum wells and solid and thin dotted curve for graphene 
layers. 



The second term in Eq. (|20[l is the temperature dependent normal density taking into account gas of phonons ("bo- 
golons") with dispersion law e{p) = fi/nigp, pi is given by Eg. p^ . 

In a 2D system, superfluidity of magnetobiexcitons appears below the Kosterlitz-Thouless transition temperature 
Tc = 7rn5(r)/(2m^), where only coupled vortices are present--.. Employing ns{T) for the superfluid component, we 
obtain an equation for the Kosterlitz-Thouless transition temperature Tc with solution 




'32 Um'^gkBT^ 

27 \ TTh^nbex 



1/3- 



1 - 1 



710 

^ c 
21/3 



(21) 



Here, is an auxiliary quantity, equal to the temperature at which the superfluid density vanishes in the mean-field 

approximation, i.e., ns{T^) = 0, = {2TTh^ncjm%/ {3C{i)))^^^ ■ The temperature = T°{B,D) may be used 
to estimate the crossover region where local superfluid density appears for magnetobiexcitons on a scale smaller or of 
the order of the mean intervortex separation in the system. The local superfluid density can manifest itself in local 
optical or transport properties. The dependence of Tc on the density of magnetoexcitons at different magnetic field 
B for superlattice consisting of quantum wells and graphene layers is represented on Fig. [31 



V. DISCUSSION 



It is shown that the low-density system of indirect magnetoexcitons in a slab of superlattice consisting of alternating 
e and h in QWs or GLs in high magnetic field occur to be instable due to the attraction of magnetoexcitons with 
opposite dipoles at large distances. Note that in spite of both QW and graphene realizations are represented by 
completely different Hamiltonians, the effective Hamiltonian in a strong magnetic field was obtained to be the same. 
Moreover, for N excitons we have reduced the number of the degrees of freedom from 2N x 2 to x 2 by integrating 
over the coordinates of the relative motion of e and h. The instability of the ground state of the system of interacting 
two-dimensional indirect magnetoexcitons in a slab of superlattice with alternating electron and hole layers of both 
QWs and GLs in high magnetic field is claimed due to the attraction between the indirect excitons with opposite 
directed dipole moments. The stable system occurs to be indirect quasi-two-dimensional magnetobiexcitons, consisting 
from indirect excitons with opposite directed dipole moments. The stability of the system is connected with the 
quadrupole-quadrupole repulsion of indirect magnetoexcitons. So at the pumping increase at low temperatures the 
excitonic line must vanish and only magnetobiexcitonic line survives. The Kosterlitz-Thouless transition to the 
superfluid state is calculated for the system of indirect magnetobiexcitons. According to Eq. (PT]) . the temperature Tc 
for the onset of superfluidity due to the Kosterlitz-Thouless transition at a fixed magnetobiexciton density decreases 
as a function of magnetic field B and interlayer separation D. This is due to the increased effective magnetic mass 
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niB of magnetoexcitons as a functions of B and D. The Tc decreases as B~^/'^ at D <C rs or as when D ^ tb- 
According to Fig. [31 the Kosterhtz-Thouless temperature is higher for the superlattice consisting of graphene layers 
than for the superlattice consisting of the quantum wells, and this difference is as stronger as the magnetic field is 
smaller. 
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